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The Milky Way Galaxy contains a large, spherical component which is believed 
to harbor a substantial amount of unseen matter. Recent observations indirectly 
suggest that as much as half of this "dark matter" may be in the form of old, very 
cool white dwarfs, the remnants of an ancient population of stars as old as the 
Galaxy itself. We conducted a survey to find faint, cool white dwarfs with large 
space velocities, indicative of their membership in the Galaxy's spherical halo 
component. The survey reveals a substantial, directly observed population of old 
white dwarfs, too faint to be seen in previous surveys. This newly discovered 
population accounts for at least 2% of the halo dark matter. It provides a natural 
explanation for the indirect observations, and represents a direct detection of 
Galactic halo dark matter. 

Dark matter in the spherical halo of the Milky Way Galaxy has been inferred because the 
gravitational field due to the known distribution of luminous matter, primarily stars, cannot 
explain the observed rotational characteristics of the Galaxy's spiral disk. A substantial 
portion of this unseen matter may be old, very cool white dwarfs {1-4)- A white dwarf is the 
extremely dense end-state in the evolution of stars with masses less than about eight times 
the mass of the Sun (M©). Once a star becomes a white dwarf, it no longer produces energy 
through nuclear fusion, and therefore cools and fades. The first four examples of ultracool 
white dwarfs, whose temperatures are below 4000 K, were identified only in the past 2 years 
(2, 3, 5-7), principally because they were too faint to detect in previous surveys. They 
have spectral energy distributions that differ dramatically from those of the hotter white 
dwarfs, consistent with white dwarf atmosphere models (8-10). The difference is due to the 
formation of H 2 molecules in white dwarf atmospheres with effective temperatures below 
~ 4500 K. Although H 2 molecules are symmetric and thus have no dipole moment, in the 
high densities of white dwarf atmospheres, collisions between the molecules are common. 
These collisions induce momentary dipole moments, which produce opacity at wavelengths 
longer than 0.6 jum (11). We conducted a survey to search for nearby, high-proper-motion 
white dwarfs that could be halo members and might exhibit this opacity. Regardless of their 
spectral appearance, nearby halo stars can be distinguished from disk stars by their space 
motion, because the two populations possess different kinematic properties. 
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Observations 

Our survey used digitized, photographic plates in the R59F and I?j passbands from the 
SuperCOSMOS Sky Survey (12-14)- R59F and Bj roughly correspond to optical wavelengths 
of 0.59 and 0.45 /im respectively. We searched for objects with proper motions, /i, between 
0.33 and 10.00 " year" 1 as faint as R59F = 19.8, using plates near the South Galactic Cap 
(SGC) at three epochs in each field. The SGC plates comprise 196 fields (15). The non 
overlapping area of each field is 25 square degrees. Hence, the area of our survey is 4900 
square degrees, or about 12% of the sky. Because of image blending problems and large 
halos of scattered light around bright stars, the actual area surveyed is 4165 square degrees, 
or 10% of the sky. This is several times larger than the areas searched in two previous 
surveys (3, 16). At the very faint end of our sample, we found many objects not included in 
the Luyten catalogs of proper motion stars (IT), primarily because Luyten's survey in the 
Southern Hemisphere did not attempt to find objects as faint as those we have found. These 
previously undetected stars are the principal members of our sample of halo white dwarfs. 

We used the technique of reduced proper motion (RPM) (18, 19) to identify white dwarfs 
by their subluminosity in comparison to main-sequence stars and their high intrinsic space 
motions. RPM is defined as Hr = R59F+ 5 log fi + 5, and is an estimate of the absolute 
magnitude based on the proper motion, \i, in arcseconds per year. We selected all 126 objects 
in the RPM diagram (Fig. |l]) having either > 22.5 or lying in the subluminous regions 
as potential halo white dwarfs. Sixty-three of these were previously cataloged (12, 20). Of 
these 63 known, high-proper-motion stars, 29 had no published spectroscopic follow-up. This 
yielded an observing list of 92 stars for follow-up, including 63 new discoveries. 

Over four nights (11 to 14 October 2000) at the Cerro Tololo Interamerican Observatory's 
V. M. Blanco 4-m Telescope, we obtained spectra in the wavelength range 0.48 to 0.98 /im for 
69 of the 92 candidates (filled symbols in Fig. |l|). We inspected the spectra within minutes of 
collecting the data, so we were able to discern by the middle of the second night which types 
of objects populated which parts of the RPM diagram (Fig. |l|). For example, as expected, 
we found only white dwarfs with Ha features on the far left of the subluminous sequence. 
Dwarf and subdwarf M-type stars are restricted to the right, which is also the bottom of 
the main sequence. We concentrated our observing resources on the cool objects, which 
occupy the lower part of the subluminous sequence. Of the 69 candidate objects observed, 
16 are M-dwarfs or M-subdwarfs, 2 are hot He white dwarfs, 13 are white dwarfs with Ha 
features, and 38 are new cool white dwarfs (Fig. §a), a few of which are probably cooler than 
WD0346+246, the prototypical ultracool white dwarf, which exhibits the H 2 opacity (2, 7, 
21, 22). 

Two of the 38 new cool white dwarfs have unusual spectra (Fig. ^o). One, LHS 1402, has 
a spectrum very similar to those of the peculiar stars LHS 3250 (5, 7) and SDSS 1337+00 
(6), but with a steeper slope toward longer wavelengths, suggesting a cooler temperature. 
The other object, WD2356— 209, possesses a bizarre spectrum, incomparable to any other 
known object. We reanalyzed the data on this object several times and found no evidence 
for residual instrumental effects. 
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White Dwarfs in the Galactic Halo 

The 38 cool white dwarfs in our sample show no spectral lines (e.g., Fig. |2]a), preventing 
determination of their radial velocities. Nevertheless, we can estimate their space motions 
from the tangential velocity, t> tan = fid, where \i is the proper motion and d is the distance 
between each star and Earth. To estimate distances, we used a photometric parallax relation 
derived from a linear least-squares fit to the cool white dwarf sample of Bergeron, Ruiz, and 
Leggett (23), supplemented by the ultracool WD0346+246, yielding the absolute magnitude 
in the Bj filter, M-q } = 12.73 + 2.58(5j — R59F). Our calculations of the spectra of very 
cool white dwarfs also show a linear relation between Mbj and Bj — R59F. The color turnoff 
in Sj — R59F due to collision-induced absorption by H2 molecules in the pure hydrogen 
atmospheres becomes apparent in Bj — R59F only at effective temperatures below 2500 K. 
The distances derived by this method are listed in Table 1. The scatter of the data around 
the least-squares fit results in a 20% uncertainty in our distance estimates. For LHS 147 and 
LHS 542, the only stars in Table 1 with measured distances (16, 24), our estimates and the 
measured distances agree within 1 and 20%, respectively. 

In order to use the velocities of the stars as an indicator of membership in the Galactic 
Halo, we transformed t> ta n the components U, V, and W, in galactic coordinates. U 
is radial away from the Galactic center, V is in the direction of rotation, whereas W is 
perpendicular to the Galactic disk. The calculations used to derive these velocities take 
into account the deviation of the velocity of the Sun with respect to the average velocity 
of nearby stars in the Galactic disk. Because our survey was centered around the SGC, 
our observations are most sensitive to the U and V velocities (Fig. Therefore, the most 
reasonable assumption is that W is zero for all of our stars. Zero is also the average value of 
W for all stars in the Milky Way. The point where U, V, and W are all equal to zero defines 
the so-called "local standard of rest," a frame of reference which is rotating in the disk about 
the Galactic center at the average rotation velocity of the disk near the Sun. Halo objects, 
because they have an isotropic velocity distribution with little or no overall rotation, should 
lag behind the local standard of rest with a distribution centered near V = —220 km s^ 1 , 
because the local standard of rest orbits the Galaxy at 220 km s _1 (25). 

Our survey contains 38 stars that unambiguously possess halo kinematics (Fig. |). We 
selected a sample of halo stars by accepting all objects whose velocities exceed the 2a = 94 
km s _1 velocity dispersion for the old disk stars (25). Our velocity cut excludes some halo 
stars with lower velocities, but permits only minimal contamination of the sample by disk 
stars. By choosing stars above 2a = 94 km s -1 , we excluded 95% of the disk stars. 

Furthermore, if a portion of our sample contains binary stars [the fraction of white dwarfs 
that are in white dwarf- white dwarf pairs is between 2 and 20% (26)}, we expect that the 
velocity distribution of halo stars (Fig. ||) will be biased toward the point (V,U) = (0,0). 
This is because we assumed that all of the stars are solitary in our estimate of distance. 
For those that are actually binary, we underestimated the intrinsic luminosity by, at most, 
a factor of 2 and the distance by a factor of \/2. Therefore the values of U and V are 
underestimated by the same factor of y2. This bias is apparent in our sample (Fig. |3|), 
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implying that while we may have included a few disk stars in our sample, we may also have 
excluded a few binary halo stars. 

Because we searched for stars having fi > 0.33" year" 1 there is a lack of stars with very 
low space motions. Furthermore, because of large differences in the epochs of the plates, an 
object with a proper motion above 3" year -1 will have a displacement of at least 30" between 
the plates. This compromises the detection of objects with extremely high space motions: 
even though we searched for objects having annual motions of up to 10", our survey probably 
had less than a 10% chance of finding very faint stars moving by 3" year -1 or more. These 
biases against stars with very low and very high space motions result in the dearth of points 
within 50 km s -1 of the point (U, V) = (0, 0) and the paucity of points at the extreme left of 
the UV velocity diagram (Fig. |^). Therefore, we underestimated the numbers of halo white 
dwarfs, and we measured a lower limit to the space density of these objects. 

The subset of objects with velocities outside the 2a velocity dispersion of the old disk 
stars (Fig. 3) contains 24 cool white dwarfs and 14 white dwarfs with Ha features (Table 1). 
An initial estimate of the space density of these objects is obtained using only the 12 coolest 
white dwarfs, whose average absolute magnitude, M R5 g F , is 16. With the survey depth of 
19.8 and the survey area, a rough volume is computed yielding a space density of 2 x 10 -4 
pc -3 , or 1.2 x 10 -4 M pc -3 . 

To obtain a more precise estimate of the space density of all 38 of the halo white dwarfs 
in our survey, we used the 1/V max technique (27), which determines the maximum volume, 
Vmax, in which the survey could have found each star, given the brightness detection limit. 
The sum of 1/V max for all of the stars in the sample is the number density. Because our 
survey was limited to stars brighter than i?59Fi im = 19.8, we have computed V maiX using the 
absolute R59F magnitudes, M R59F , implied by the distances inferred above. 

Knax = (4/3)7^0/(470 PC 3 (1) 

Here, Q is the surveyed area in steradians, and <i max , the maximum distance in parsecs (pc) 
that determines V mauX , is the minimum of the two following relations. 

log rf max (pc) = 0.2(i?59F lim - M R59F ) + 1 (2) 
dmax (pc) = d/V0.33 (3) 
The number density derived in this manner is 1.8 x 10 -4 pc -3 . Using 0.6 M© as the 
average white dwarf mass in the halo (28), the local space density of halo white dwarfs is 
1.1 x 10 -4 M & pc -3 . However, our survey is not complete to the R59F = 19.8 level. A 
measure of the uniformity of our sample is given by the average value of V/V maK , where V is 
the volume of space given by the distance to the given star and the survey area. In the case 
of this survey, the average of V/V max is 0.46. For a sample distributed uniformly in space, 
this value should be 0.5. An average of V/V max = 0.5 is found for R59F = 19.7, indicating 
that all stars brighter than 19.7, with the proper motion constraint mentioned above, were 
found in the survey. Recalculating the space density using _R59Fi im = 19.7 yields 1.3 x 10 -4 
M pc -3 . 

The density of white dwarfs in the halo predicted from the subdwarf star counts and a 
standard initial mass function is 1.3xl0 -5 M Q pc -3 (29), which is ten times smaller than the 
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value we calculate. This means that star formation in the early Galaxy must have favored 
higher-mass stars that would evolve into the white dwarfs we are detecting now. At the same 
time, this early star formation must have under produced the low-mass subdwarfs relative 
to the more recent star formation processes observed in the disk (4). 

The estimated density of halo dark matter near the Sun is approximately 8xl0~ 3 M Q 
pc -3 (30). Thus, the population we have detected accounts for about 2% of the local dark 
matter. We treat this number as a lower limit for the reasons discussed above and because 
it appears that we have not probed the entire range of intrinsic luminosities of these objects. 
Indeed, the number density as a function of M R59F suggests that we are only detecting the 
rising power law of the luminosity function and have not seen any indication of the expected 
turnover due to the finite life time of the halo white dwarf population (31). This means that 
there may be an undetected, larger population of even fainter and cooler white dwarfs in the 
halo. 

Discussion 

Microlensing experiments have indirectly revealed a population of massive, compact ob- 
jects in the line of sight to the Large Magellanic Cloud. The Massive Compact Halo Object 
project (MACHO) (1) estimates that these compact objects contribute between 8 and 50% of 
the local halo dark matter at the 95% confidence level, and the inferred lens mass suggests 
that they are white dwarfs. The similar Experience pour la Recherche d'Objets Sombres 
project (EROS) places a strong upper limit of 35% for such objects (32). Because the halo 
white dwarf space density that we derive from our survey represents a lower limit and is con- 
sistent with the two measurements mentioned above, the population of white dwarfs with 
halo kinematics that we have discovered provides a natural explanation for the microlensing 
results. 

A number of the white dwarfs revealed in this survey may be cooler than any other white 
dwarf previously known. There are at least three objects (Table 1) that are comparable 
in temperature to WD0346+246 (7), representing a doubling of the number of ultracool 
white dwarfs known. The nature of cool white dwarfs can be discussed on the basis of their 
positions in color-color diagrams. Although a combination of optical and infrared colors is 
necessary at these low effective temperatures (7, 23), for now, we must contend with the 
limited spectral coverage of photographic plates for this sample. The comparison of Bj, 
R59F, and In (a bandpass centered on 0.8 /jm) colors is a poor diagnostic of atmospheric 
composition for effective temperatures, T e fj, greater than 3500 K (Fig. f|). Given the large 
photometric error bars, there is little that distinguishes the colors of these stars from a 
sample of cool white dwarfs with disk kinematics (33). 

Stars with Bj—R59F ^1.2 are relatively warm, with T e fj ^ 5000 K (except for the strange 
stars LHS 3250, SDSS 1337+00 and LHS 1402). They are also the most luminous and most 
distant stars in the sample (Table 1). Although it is not possible to guess the atmospheric 
composition of any of these stars from the present data, we point out that the sample may 
contain a few unusual objects, in addition to the peculiar LHS 1402 and WD2356— 209 (Fig. 
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2b). Other inferences, however, can be made from the photographic color-color diagram (Fig. 
f|), using stars previously studied. LHS 542 has a pure He atmosphere and T e s = 4747 K (24)- 
WD0346+246 and F351— 50 are similar to each other with predominantly He atmospheres, 
a very small admixture of H, and T e g- of ~ 3750 K and ~ 3500 K, respectively (7). The 
most extreme star of the sample, and also the faintest, is WD0351— 564 with Bj — R59F = 
1.98. The stars with B 3 - R59F > 1.6 (WD0351-564, WD0205-053 and WD0345-362) 
may have either pure He atmospheres with T e g > 4000 K or, more interestingly, be other 
examples of very cool white dwarfs with He atmospheres polluted by a small amount of H, 
like WD0346+246. The stars most likely to have pure H atmospheres residing below the 
color turnoff due to strong H 2 collision-induced opacity are WD0135— 546, WD0340— 330, 
and, perhaps, WD0351— 564. 

We still do not understand the nature of the three strange objects LHS 3250, LHS 1402, 
and SDSS 1337+00 (5-7), all of which lie at R59F — In colors near —0.5 (well outside of 
Fig. |4j). The kinematics of these stars, with the exception of LHS 1402, suggest that they 
may be disk objects. Although they are certainly cool white dwarfs, we know neither their 
effective temperatures nor their atmospheric compositions. They may represent another, 
perhaps unusual, stage in the spectral evolution of cool white dwarfs. 

The direct detection of a significant white dwarf halo population opens promising avenues 
of investigation. Besides the opportunity to study individual white dwarfs that are older and 
cooler than any known so far, dedicated studies of this population of stellar remnants will 
bear directly on the ancient history of the Galaxy. For instance, a determination of the 
luminosity function of halo white dwarfs and a characterization of the low luminosity cutoff 
will reveal the age of the first generation of halo stars, lead to a determination of its initial 
mass function, and provide strong constraints on formation scenarios of the Galaxy. The 
interpretation of observations of young galaxies at high redshifts will benefit from a better 
understanding of the early history of our Galaxy. Finally, the complete characterization of 
this new stellar population will reveal whether it is responsible for the microlensing. 
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Table 1. Candidate Halo White Dwarfs. 
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53 


WD0214-419 


02:14:14.887 


-41:51:09.04 


0.334 ±0.026 


106.7 


20.63 


1.17 


0.48 


150 


95 


WD0044-284 


00:44:02.188 


-28:24:11.15 


0.348 ± 0.063 


191.1 


20.83 


1.48 


0.55 


119 


72 


WD2214-390 C 


22:14:34.727 


-38:59:07.05 


1.056 ±0.013 


110.1 


16.57 


0.75 


0.40 


120 


24 


WD2324-595 C 


23:24:10.165 


-59:28:07.95 


0.581 ± 0.016 


167.1 


17.13 


0.23 


0.05 


159 


58 


LP 588-37 c 


01:42:20.770 


-01:23:51.38 


0.354 ±0.010 


162.8 


18.62 


0.19 


0.0 b 


202 


120 


WD0345-362 


03:45:32.711 


-36:11:03.60 


0.589 ±0.081 


167.8 


21.12 


1.81 


0.60 


155 


55 


WD0045-061 


00:45:06.325 


-06:08:19.65 


0.699 ±0.017 


171.2 


18.82 


1.12 


0.66 


145 


44 


WD0225-326 


02:25:28.681 


-32:37:53.92 


0.338 ± 0.033 


59.7 


18.94 


0.58 


0.13 


140 


88 


WD2348-548 


23:48:46.904 


-54:45:46.21 


0.388 ± 0.033 


105.5 


19.76 


1.14 


0.51 


121 


66 


WD01 17-268 


01:17:51.649 


-26:48:51.21 


0.468 ± 0.032 


82.2 


19.91 


1.34 


0.58 


122 


55 


LP 651-74 c 


03:07:14.113 


-07:14:59.12 


0.477 ± 0.036 


203.8 


17.56 


0.94 


0.4 b 


68 


30 


WD0135-546 


01:35:38.677 


-54:35:27.78 


0.675 ± 0.018 


80.9 


19.52 


1.47 


0.43 


127 


39 


WD0100-567 C 


01:00:43.076 


-56:46:36.61 


0.414 ± 0.043 


45.1 


17.85 


0.63 


0A b 


98 


50 



"This object is extremely cool but has a spectrum similar to the peculiar LHS 3250 (7). 
6 These R59F — 1^ colors are calculated from spectrophotometry. 
c These stars have Ha features. All others have featureless spectra. 



Table 1: Candidate Halo White Dwarfs. WD names are new discoveries; otherwise, Luyten 
LHS or LP names are given; J0014-3937 was discovered by Scholz et al. (20), who classify 
this object as "cool no Ha", and F351-50 was discovered by Ibata et al. (3). The tangential 
velocity, t> ta n(est) in km s -1 , is an estimate based on photometric parallax. Distances are 
estimates based on a color- magnitude relation (see text). Star positions are given in right 
ascension (RA) and declination (Dec.) at the equinox and epoch of J2000. The proper 
motion, /i, and position angle, PA, are given in " year" 1 and degrees, respectively. The 
photometric measurements in Bj, R59F, and J N are given in magnitudes, with uncertainties 
of ±0.1 mag. on average. The distance is in parsecs from the Sun. 
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Figure 1: Reduced proper motion diagram. All low-luminosity, proper-motion objects in 
our survey are shown. The main sequence extends from the top-center downward. The 
subluminous sequence to the left and below is the white dwarf sequence. Blue, filled stars 
indicate white dwarfs confirmed with our spectroscopy or previously known to be white 
dwarfs. Open stars are objects we have not yet observed spectroscopically but presume to 
be white dwarfs on the basis of their location in this diagram. Filled, red triangles are M- 
dwarfs and filled, green circles are subdwarfs for which we obtained spectra. Open triangles 
and open circles are suspected M-type dwarf or sub-dwarf stars, respectively, which we did 
not observe. The dashed line indicates the value of Hr below which the Luyten catalogs (17) 
contain very few objects. 



10 




A (/Oim) A (/u,m) 

Figure 2: Selected spectra of new cool white dwarfs. WD0346+246 is included for compari- 
son, (a) contains five spectra representative of the set of 24 cool halo white dwarfs (without 
Ha features) such that the full range of spectral shapes in this set is presented. The spec- 
tra, normalized and shifted vertically to facilitate comparison, are ordered according to the 
sequence in Fig. |], from top down in order of increasing Bj— R59F. This sequence roughly 
corresponds to decreasing temperature, (b) shows the odd spectrum of WD2356— 209, which 
has no analogs, the spectrum of LHS 1402, which seems to be a cooler analog of LHS 3250 
and SDSS 1337+00, and a 4000 K blackbody spectrum (BB) at the top. Object names used 
in Table 1 are indicated. None of the fine features are real. 
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Figure 3: Velocities of white dwarfs in the survey. The terms U and V are the Galactic radial 
and rotational velocities, respectively, of the white dwarfs in our survey Objects that appear 
in Table 1 are surrounded by boxes. The dashed ellipses indicate the velocity dispersions 
of the Galaxy's old disk (right) and halo (left). The solid ellipses are the 2a dispersions for 
these populations. The ellipses representing the old disk stars, the highest velocity members 
of the Galaxy's disk, are centered at (V,U) = (—35,0), because this population of stars 
rotates slightly slower than the local standard of rest. 
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Figure 4: Photographic color-color diagram. The cooler stars listed in Table 1 are shown as 
solid dots, and a typical error bar is indicated. The synthetic colors of pure He atmospheres 
(open squares) (9) form a diagonal sequence from T c s = 8000 K (lower left) to 4000 K, in 
steps of 1000 K. For pure H atmospheres (open triangles) (9) the sequence starts at 7000 K 
(lower left) and turns over to bluer R59F — I N below 3500 K due to the onset of H 2 collision- 
induced opacity. The dotted line connecting open circles shows the effect of mixed H/Hc 
composition on a T cS = 3500 K model, with the number ratio of He to H ranging from to 
10 7 (7). All models have logg(cm s~ 2 )= 8. As a basis for comparison, three stars whose 
spectra have been studied in greater detail (3, 7) are labeled. WD0346+246 is not part of 
our survey and is shown here for comparison. 



13 



